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We report quantum dynamics calculations of the O + OH ^ H + O2 reaction on two different rep- 
resentations of the electronic ground state potential energy surface (PES) using a time-independent 
quantum formalism based on hyperspherical coordinates. Calculations show that several excited 
vibrational levels of the product O2 molecule are populated in the reaction. Rate coefficients eval- 
uated using both PESs were found to be very sensitive to the energy resolution of the reaction 
probability, especially at temperatures lower than 100 K. It is found that the rate coefficient re- 
mains largely constant in the temperature range 10 - 39 K, in agreement with the conclusions of a 
recent experimental study [Carty et al., J. Phys. Chem. A 110, 3101 (2006)]. This is in contrast 
with the time-independent quantum calculations of Xu et al. [J. Chem. Phys. 127, 024304 (2007)] 
which, using the same PES, predicted two orders of magnitude drop in the rate coefficient value 
from 39 K to 10 K. Implications of our findings to oxygen chemistry in the interstellar medium are 
discussed. 

PACS numbers: 



I. INTRODUCTION 



The reaction 



0(=^P) + 0H(2n) ^ H(2S) 



) 



is of considerable importance in atmospheric, combustion 
and interstellar chemistry. In the upper stratosphere and 
lower mesosphere, it plays a major role in partitioning the 
relative abundance of OH and HO2 molecules [H, S [I]- 
The reaction also plays a key role in the night-time air- 
glow emissions from the hydroxyl radical which has also 
been recently detected in the atmosphere of Venus [1] . It 
has been identified as a key reaction in interstellar oxy- 
gen chemistry [1, d, 01 and it is considered to be the most 
important source of oxygen molecules in cold interstellar 
clouds [1, 0. As a consequence, the HO2 system has 
been the topic of numerous electronic structure calcula- 
tions of its potential energy surface p^. [Til. [T^. [Tsl. [l^. [isj 
as well as dynamics calculations to evaluate its temper- 
ature dependent rate coefScients [H, [13, [H, d) Elj 
HO, mi 111 m, [13, iS iij. The reaction has also been 
the focus of a large number of experimental measure- 
ments [13, H, m [13, [U, m, m, H il Hi- How- 
ever, there still remains significant discrepancy between 
measured and computed values of its rate coefficients for 
temperatures below 200 K, which is the most important 
region for astrophysical and upper atmospheric applica- 
tions. 

Recently, Carty et al. [SG*] reported an experimental 
measurement of rate coefficients for the reaction in the 
temperature range 39 K to 142 K and observed no vari- 
ation in the rate coefficients with temperature in this 
regime. Based on their findings they concluded that the 
rate coefficient would largely remain constant between 
39 K to 10 K, temperatures typical of cold interstellar 



clouds. In contrast to the experimental results of Carty et 
al., in a recent theoretical work using a time- independent 
quantum mechanical (TIQM) method and a J-shifting 
approximation, Xu et al. [2J] reported a rate coefficient 
that precipitously dropped by about two orders of mag- 
nitude between 39 K and 10 K. Using a time-dependent 
wave packet (TDWP) method that does not employ the 
J-shifting approximation, Lin et al. [25| reported a rate 
coefficient at 10 K that is about one order of magnitude 
smaller than that obtained by Carty et al. at 39 K. In 
a subsequent study, Quan et al. [371] adopted these rate 
coefficients for modeling molecular oxygen chemistry in 
the interstellar medium. 

In this paper, using an accurate TIQM method and two 
different representations of the electronic ground state of 
the HO2 system, we show that the low temperature rate 
coefficient of the O -|- OH reaction is very sensitive to the 
dense resonance structures in the energy dependence of 
its reaction probability. By using a very fine grid of colli- 
sion energies in the low and ultralow regime to compute 
the reaction probabilities, we obtain nearly temperature 
independent values of the reaction rate coefficients in the 
range 10 - 39 K, in agreement with the conclusions of 
Carty et al. [s^ . Our results differ from those of Xu et 
al. [2JI which show rapid decrease below 40 K. However, 
our results do merge with those of Xu et al. at tempera- 
tures above 300 K where the low energy regime does not 
make a significant contribution. We believe that the rate 
coefficient reported here would provide a more accurate 
value for modeling oxygen chemistry in the interstellar 
medium. 

A brief description of the computational method is pro- 
vided in section II followed by results and discussion in 
section HI. Conclusions are presented in section IV. 
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II. COMPUTATIONAL METHOD 



The calculations have been performed using the adia- 
batically adjusting principal-axis hyperspherical (APH) 
approach of Pack and Parker [s^. In the present work, 
we employed two representations of the electronic ground 
state (1 ^A") of the HO2 system. We used the XXZLG 
PES calculated by Xu, Xie, Zhang, Lin, and Guo [l3. ITp 
which was also employed in the study of Xu et al. [2J|. 
For the other surface, we used the diatomics-in-molecule 
(DIM) PES developed by Kendrick and Pack [l3|, here- 
after referred to as the DIMKP PES. For each value of the 
total angular momentum quantum number, J, and each 
value of the hyperspherical radius, p, the wavefunction is 
expanded onto a basis set of adiabatic functions, which 
are eigenfunctions of a triatomic hyperangular Hamil- 
tonian. The hyperangular Hamiltonian is diagonalized 
using an Implicitly Restarted Lanczos algorithm and a 
hybrid DVR/FBR primitive basis set [3^. The time- 
independent Schrodingcr equation yields a set of dif- 
ferential close-coupling equations in p, which is solved 
using the log-derivative matrix propagation method of 
Johnson [i^l • The log-derivative matrix is propagated to 
a matching distance where asymptotic boundary condi- 
tions are applied to evaluate the reactance matrix, K'^ , 
and the scattering matrix, S"' . The square elements of 
the S"^ matrix provide the state-to-state transition prob- 
abilities, P"^ . The matching distance pm and the number 
of adiabatic functions n used in the basis set are deter- 
mined by optimization and extensive convergence stud- 
ies. The convergence of the J = reaction probability 
Pv=o^j=o(^c) for the O + OH(t; = 0, j = 0) reaction with 
respect to pm and n is summarized in Table |T] for the 
DIMKP PES and Table lU for the XXZLG PES. Based 
on these convergence studies, we have used pm — 26.8 oq 
and n = 393 in the final production calculations. 



pm (ao) 



26.8 
393 



32.7 

393 



160.6 

393 



160.6 
380 



160.6 
300 



(eV) P„lt,^o 

0.001 0.7269 0.7295 0.7284 0.7284 0.7282 

0.01 0.3498 0.3492 0.3490 0.3490 0.3484 

0.1 0.2956 0.2956 0.2954 0.2954 0.2956 

1 0.1831 0.1841 0.1825 0.1835 0.2653 



pm (ao) 
n 



26.8 
393 



32.7 
393 



160.6 
393 



160.6 160.6 
380 300 



(eV) P:Lt!=o 
0.001 0.2153 0.2007 0.2081 0.2081 0.2076 
0.01 0.2563 0.2570 0.2577 0.2577 0.2568 
0.1 0.1816 0.1807 0.1795 0.1795 0.1792 

1 0.0846 0.0846 0.0841 0.0846 0.1355 



TABLE 11: Same as Table m but for the XXZLG PES. 



III. RESULTS AND DISCUSSION 



Probabilties 



In Fig. [T] we show the J = reaction probabilities 
as a function of the collision energy computed using the 
DIMKP PES (upper panel) and the XXZLG PES (lower 
panel) . The reaction probabilities calculated by Xu et al. 



DIMKP PES 




TABLE I: J = reaction probability of O -I- OH(ii = 0, j = 
0) R + O2 for the DIMKP PES for different values of 
the matching distance pm and the number of hyperspherical 
channels n at different collision energies Ec. 



FIG. 1: (Color online) J = reaction probabilities of O 
4- OH(w = 0,j = 0) H -I- O2 reaction on the DIMKP 
PES (upper panel) and the XXZLG PES (lower panel). As 
discussed in results of Xu et al. ^ on the XXZLG PES 
have been shifted back by 0.00409 eV to make a one-to-one 
comparison with our results. 
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O + OH(v=0J=0) ^ H + Ojiv,) ; DIMKP PES 




(eV) (eV) 



FIG. 2; Final vibrational distribution of the J — reaction 
probability of O + OH(w = 0, j = 0) H + 02(u/) for the 
DIMKP PES. 



using the XXZLG PES and a TIQM method are also in- 
cluded in the lower panel. In both figures, the inset shows 
results for Ec < 0.03 eV to illustrate the energy resolu- 
tion required for resolving the resonance features in the 
reaction probability. Our calculations include reaction 
probabilities in the energy range Ec — [10^^ — 0.8] eV. 
In the energy range Ec = [10~^ — 10"*] eV where no res- 
onance features are present we included 30 energies in a 
logarithmic scale. Above 0.0001 eV, we used the follow- 
ing energy grids: Ec = [0.0001 - 0.0600; 0.0001] eV and 
Ec = [0.060-0.800; 0.001] eV using linear intervals where 
the last number in the brackets indicates the energy spac- 
ing. The reaction probabilities obtained using the two 
PESs are not identical and they illustrate the sensitiv- 
ity of results to details of the interaction potential. The 
global mean is of about 0.3 for the DIMKP PES while it 
is about 0.2 for the XXZLG PES. The overall trend of the 
reaction probabilities for the XXZLG PES is a decrease 
with increase in the collision energy while it oscillates 
around a value of about 0.3 for the DIMKP PES. The 
large number of resonances seen in the reaction probabil- 
ities comes from quasi-bound states of the HO2 complex. 
We note that our results on the XXZLG PES are in ex- 
cellent agreement with those of Xu et al. (red curves) at 



O + OH(v=0,j=0) ^ H + Ojlv,) ; XXZLG PES 




E^ (eV) (eV) 
FIG. 3: Same as Fig. [3 but for the XXZLG PES. 



high and intermediate collision energies. However, some 
discrepancies appear for low energies Ec < 0.012 eV as 
seen in the inset of the figure. As discussed in the next 
section, these discrepancies at the lower energies can lead 
to significantly different rate coefficients at low tempera- 
tures. 

Since the O + OH(w = 0, j = 0) ^ H + OaCw/) re- 
action is exoergic, the reaction populates a number of 
excited vibrational levels of the O2 molecule. The rela- 
tive kinetic energy of the ensuing products depends on 
the vibrational population of the O2 molecule. Since 
non-equilibirium kinetics is an important issue in upper 
atmospheric and astrophysical environments nascent vi- 
brational populations of the O2 molecules resulting from 
the O + OH reaction will be important in modeling hy- 
droxyl and oxygen chemistry in these environments. The 
computed vibrational distributions are shown in Fig. [2] 
for the DIMKP PES and in Fig.[3]for the XXZLG PES. 
The results for the DIMKP PES show that for coUision 
energies Ec < 0.4 eV, O2 molecules are formed preferen- 
tially in low- lying vibrational levels Vf — — 3. These 
levels are open for the O -I- OH reaction, even in the 
limit of vanishing collision energies. For collision ener- 
gies Ec > 0.4 eV, formation of O2 molecules in Vf = 4, 5 
also competes with the lower vibrational levels. The re- 
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suits for the XXZLG PES in Fig. [3] show similar trends. 
Vibrational levels vj — — 3 are the most probable for 
Ec < 0.4 eV, while Vf — 4,5 are equally probable for 
Ec > 0.4 eV. For Vf = 0-3 and Ec < 0.4 eV, the reaction 
probabilities on the XXZLG PES are larger than that of 
the DIMKP PES. The opposite is true for Ec > 0.4 eV. 
Overall, the reaction probabilities for a given Vf are less 
sensitive to the collision energy for the DIMKP PES than 
that of the XXZLG PES, consistent with the result for 
the total reaction probabilities shown in Fig. [T] 



B. Rate coefficients 

Accurate determination of the rate coefficients would 
require calculations of the reaction probabilities for all 
contributing values of J. Computational expense esca- 
lates quickly with J unless some angular momentum de- 
coupling approximations are used. Like Xu et al., we use 
the J-shifting approximation (4^ to compute the initial 
state-selected rate coefficients for the reaction. Within 
the J-shifting approximation, the rate coefficient is given 
by the expression: 



X P::'=\Ec) e-^^/C^--^) dEc (1) 



where fc^ is the Boltzmann constant, P^'j ^ is the reac- 
tion probability and E'^YiUt the height of the effective 
barrier for a given partial wave J in the entrance chan- 
nel. The barrier height is evaluated from the effective 
potential, V|,g, for a given partial wave: 



*^ err 



h^J{J + i) 

2/i(i?0-OH)2 



Knin(-Ro-OH) 



(2) 



where Vinin(-Ro-OH) is the minimum energy path of the 
PES as a function of i?o-OH and /i is the O— OH reduced 
mass. The minimum energy paths for both PESs are 
shown in Fig. g] In Eq. (P), Qr = Qtrans x Qoi is the 
reactant partition function. For the translational par- 
tition function we used the standard formula, Qtrans = 

( ) ■ For the electronic partition function we used 
the expression given by Graff and Wagner [l8| : 

^ _ (5 + 3e-228/T^ + e-326/T) (2 + 2e-205/T) 



1. Sensitivity of rate coeffieients to energy resolution of 
reaction probabilities 

Calculation of rate coefficients involves integration of 
the Boltzmann distribution weighted by the reaction 
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FIG. 4: Minimum energy path of the O + OH(u = 0, j = 0) 
reaction as a function of iio-OH for the DIMKP and XXZLG 
PESs. 



probability over the collision energy. The integration has 
to be performed carefully at low temperatures where the 
Boltzmann distribution can extend to low energies. This 
is especially important for capture-type reactions as the 
present system for which the rate coefficients are gener- 
ally large at low temperatures. The typical temperature 
in the cold interstellar clouds is about 10 K. At this tem- 
perature, one needs to include collision energies lower 
than 10~^ eV (w 1 K) in accurately calculating the rate 
coefficients. It is known from Bethe-Wigner laws [45l.l46j 
that rate coefficients of exothermic reactions are finite in 
the limit of zero temperatures [47| . The limiting value 
can be large for tunnelirig dominated reactions such as 
the F -I- H2 system [U I3 ^-s well as barrierless reac- 
tions such as Li + Li2(t') collisions [i^. Furthermore, 
numerous triatomic resonances can appear in the prob- 
abilities at collision energies near the reaction threshold 
which can affect the value of the rate coefficient at low 
temperatures. 

We have paid careful attention to the convergence 
of the reaction probabilities and rate coefficients with 
the resolution of collision energy grid at very low ener- 
gies. Four sets of collision energies have been used to 
check the convergence of the rate coefficients. Fig. [3] 
shows the low-energy portion of the reaction probabil- 
ities employed in the four sets for the two PESs. Us- 
ing the same notation as before and with the energies 
in eV, the first set (set 1, bold red curve) corresponds to 
Ec = [0.001-0.800; 0.001]. The second set (set 2, dashed 
green curve) includes Ec = [0.0001 - 0.0600; 0.0005] and 
Ec = [0.060 - 0.800; 0.001]. The third set (set 3, thin 
blue curve) uses a finer energy resolution in the low 
energy regime: Ec = [0.0001 - 0.0600; 0.0001]; = 
[0.060-0.800; 0.001]. The fourth set (set 4, dashed black 
curve) is composed of set 3 plus the ultralow energy 
regime, Ec = [10^^ — 10^"*]. There is no resonance fea- 
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10' 

EJeV) 



FIG. 5: (Color online) Convergence of the reaction probabil- 
ities as functions of the collision energy for the DIMKP PES 
(upper panel) and for the XXZLG PES (lower panel). See 
text for the definition of the sets employed. 



tures in the ultralow energy regime where the probabili- 
ties vary as the square root of the energy, in accordance 
with the Bethe-Wigner laws [i^, . 

As Fig. [5] shows the ultralow regime is more prominent 
for the DIMKP PES due to the different description of 
the long-range interaction potential. It is seen that the 
resonances are not fully resolved in calculations with set 
1 and set 2 for both PESs. This is especially the case for 
the XXZLG PES in the energy range 0.001 - 0.005 eV. 
Thus, the third multiplicative term in Eq. ([T]) is more 
accurately evaluated using set 3 and set 4. 

Figure [6] shows the convergence of the rate coefficients 
with the energy resolution for the four collision energy 
sets for the DIMKP PES (upper panel) and the XXZLG 
PES (lower panel). The figure clearly illustrates that the 
rate coefficient at temperatures below 100 K are very 
sensitive to the resolution of the energy grid used. The 
convergence improves with increase in energy resolution 
of the reaction probabilities and also when the ultralow 
energy regime is included. The contribution of the ul- 
tralow energy regime is more important for the DIMKP 
PES as evident from the corresponding reaction proba- 
bilities shown in Fig. \5\ Figure [S] shows that careful at- 
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FIG. 6: (Color online) Same as Fig. [5] but for the rate coef- 
ficients as functions of the temperature. 



tention must be devoted to the resolution of the energy 
grid and/or to the inclusion of the very low energy regime 
in order to accurately calculate the low temperature rate 
coefficients of the O -t- OH reaction. This applies regard- 
less of whether or not a J-shifting approximation is used 
in the calculation of the rate coefficients. 



2. Comparison with experiments and other theoretical 
results 

The converged rate coefficients (results from set 4 
above) for the v — 0,j = initial state obtained from 
the DIMKP and the XXZLG PESs are shown in Fig. [7] 
as functions of the temperature along with experimental 
data from several groups [29', '55, '36*1 as well as the rec- 
ommended values by NASA [43.] and lUPAC [44|]. The 
theoretical results of Xu et al. [l^l obtained using the 
XXZLG PES and the results of Harding et al. ^ ob- 
tained using the PES of Troe and Ushakov [l^ are also 
included for comparison. Table Hill lists rate coefficients 
at selected temperatures from different theoretical and 
experimental studies. 

Fig. [7] shows that the XXZLG PES yields re- 
sults in somewhat better agreement with the exper- 
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FIG. 7: Rate coefficients of O + OH(i; 
O2 as functions of the temperature. 



0,j = 0) ^ H + 



iments than the DIMKP PES. At T 298 K, 

the rate coefficient calculated using the XXZLG 
PES is 3.01x10^^^ cm'^ molecule"^ s~^ compared to 
2.45x10^^^ cm'^ molecule"^ s~^ obtained using the 
DIMKP PES. The result on the XXZLG PES is within 
the reported error bars of the NASA panel recommended 



value of 3.3±0.7xl0' 



molecule 



at 298 K. 



The computed results on the DIMKP and XXZLG PESs 
lie within the quoted error bars of Carty et al. at 
T = 142 K. At T = 39 K, the result of the XXZLG is 
slightly above the experimental result, while it is higher 
for the DIMKP PES. 

The temperature dependence of the rate coefficients 
predicted by the two PESs is quite similar. They predict 
rate coefficients within about 30% for 10 < T < 100 K 
with the DIMKP PES yielding higher values. This can be 
explained by the energy dependence of the reaction prob- 
abilities shown in the insets of Fig. [1] The reaction prob- 
abilities are somewhat higher for the DIMKP PES than 
for the XXZLG PES between - 0.01 eV (« - 100 K), 
leading to an increase in the third multiplicative term in 
Eq. ^ . This can also be explained by the topology of the 
PESs. As shown in Fig. 31 the minimum energy path, 
Vmin(i?o-OH), of the DIMKP PES features a small "reef 
at i?o-OH — 5.0 gq. It is located at i?o-OH = 5.4 gq for 
the XXZLG PES. For the latter PES the reef is about 
0.08 eV higher than that of the DIMKP PES. The loca- 
tion and height of the reef is sensitive to the electronic 
structure method employed. For certain values of J, it 
becomes an effective barrier. The smaller reef for the 
DIMKP PES leads to smaller effective barriers especially 
for higher J values. This enhances the second multiplica- 
tive term in Eq. H]) and leads to a higher rate coefficient 
for the DIMKP PES. For T >100 K, the two PESs pre- 
dict rate coefficient within 20% with the DIMKP PES 
yielding smaller values. This is attributed to the smaller 
overall reaction probabilities of the DIMKP PES com- 
pared to the XXZLG PES for Ec > 0.01 eV (w 100 K), 
as seen in the insets of Fig. [T] 



Our computed rate coefficients do not show a signifi- 
cant decrease between 142 K and 39 K, consistent with 
the experimental results of Carty et al. [IBl . Though no 
experimental data are available for temperatures below 
39 K, our results on both PESs do not predict a dramatic 
decrease between 39 K and 10 K. At 10 K we obtain 
a rate coefficient of 3.91x10^^^ cm'^ molecule"^ s^^ on 
the XXZLG PES and 5.25x10"" cm^ molecule"! s^^ 
on the DIMKP PES. Our results are about a factor of 70 
larger than the value of 5.41 xlO"!'^ cm^ molecule"^ 
reported by Xu et al. [Hi using the XXZLG PES and 
the same J-shifting approximation. The discrepancy is 
attributed to the sparse energy grid, small differences in 
the reaction probabilities for Ec < 0.012 eV, and the arti- 
ficial energy shift in the calculations of Xu et al. "5^1 . 
Among these, the energy shift is the main source of the 
discrepancy. Based on our results, we believe that the 
rate coefficients calculated by Xu et al. may not be 
appropriate for modeling oxygen chemistry in the inter- 
stellar medium [33] ■ While an accurate calculation of 
the rate coefficient would require the inclusion of many 
higher angular momentum quantum numbers and non- 
adiabatic couplings, the present study shows that special 
attention must be given to the energy grid in the cal- 
culation of low temperature rate coefficients for capture 
reactions. 



reference 



10 K 39 K 142 K 298 K 



this work, DIMKP 5.25 

this work, XXZLG 3.91 

Xu et al. [24], XXZLG 0.0541 

Lin et al. [25], XXZLG 0.784 
Carty et al. [3^ 
NASA [43| 
lUPAC [44] 



5.71 
4.66 



3.90 
4.30 



2.45 
3.01 



3.5±1.0 3.5±1.0 



3.3±0.7 
3.5±0.4 



TABLE III: Rate coefficients of O -I- OH ^ H -I- O2 in units 
of 10"!! cm'^ molecule"! s"! for different temperatures. 



IV. CONCLUSION 

In conclusion, we have performed quantum dynam- 
ics of the O -I- OH — > H -|- O2 reaction over a wide 
range of collision energies on two recent PESs using a 
time-independent quantum formalism based on hyper- 
spherical coordinates. We report total and product vi- 
brational state-selected reaction probabilities, and initial 
state-selected rate coefficients for the reaction. Special 
attention has been devoted to the convergence of the rate 
coefficients with respect to the energy resolution of reac- 
tion probabilities and the inclusion of the ultralow energy 
regime. The computed rate coefficients are in reasonable 
agreement with experimental measurements and existing 
recommended values, despite using a J-shifting approx- 
imation for the evaluation of the rate coefficients. Our 
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calculations show that omission of the low energy regime 
or insufficient energy resolution of the reaction probabil- 
ities can lead to significant errors in the computed rate 
coefficients. The rate coefficient at 10 K differ by a factor 
of 2-3 between calculations using a finer energy grid and 
a sparse energy grid. Therefore, we expect that a careful 
choice of the energy grid will also be important for an 
accurate evaluation of the rate coefficients without the 
J-shifting approximation. Overall, results on the XX- 
ZLG PES are in bette r ag reement with the experimental 
results of Carty et al. ^ compared to the DIMKP PES 
in the temperature range of 39—142 K. However, only 
two data points are available from measurements in this 
temperature range and experimental error bars are also 
quite large. For both PESs, the rate coefficients decrease 
only slightly as the temperature is decreased from 39 K 
to 10 K, in agreement with the conclusions of Carty et 
al. Based on the present results we believe that a re- 
evaluation of the importance of the O -|- OH reaction on 



O2 abundance would be required for describing oxygen 
chemistry in the interstellar medium. 
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